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Statics and dynamics of a one-dimensional quantum many-body system 
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The macroscopic zero-temperature behavior of weakly- incommensurate systems in one dimension 
is described in terms of sohtons. The soliton density n obeys equations displaying several types of 
singular interface-like solutions: (i) equilibrium or moving boundary between the n = and finite n 
regions, and (ii) stationary or moving annihilation front separating solitons from antisolitons. 



Low-dimensional materials provide a testing ground 
for the theories of strongly correlated systems Q|. Re- 
cently there have been several attempts to manufacture 
truly one-dimensional systems with periodically modu- 
lated parameters. Kouwenhoven et al. Q built a one- 
dimensional wire of spinless electrons with a 200 nm pe- 
riodic potential. Upon varying the electron density with 
a gate voltage they observed suppression of conductance 
at some rational fillings. Recently Tarucha et al. ||] has 
succeeded in introducing a shorter period potential. A 
different path was undertaken by van Oudenaarden et 
al. who modelled a one-dimensional quantum liquid 
placed in a periodic potential by a very long and narrow 
array of small Josephson junctions placed in a magnetic 
field. The field determines the density of Josephson vor- 
tices in the array; these behave like quantum particles. 
A suppression of the array resistance at some rational 
fillings was observed, and interpreted in terms of a Mott 
transition. Another relevant experimental system would 
be a weakly-doped nanotube at low temperature. 

In each of these examples, the fundamental model is 
a one-dimensional quantum liquid ~ spinless fermions, or 
bosons with a short-ranged repulsive interaction - at zero 
temperature, placed in an external periodic potential . 
The possible phases are classified as follows: 

*When the interactions between particles are not 
strongly repulsive, quantum fluctuations will render the 
periodic potential irrelevant, and the system will always 
be a conductor. 

*Strongly repulsive particles will become immobilized 
even by a weak commensurate potential, and the result- 
ing state has long-ranged order, so that we may idealize 
this structure as being a crystal of particles in registry 
with the periodic potential. This state is called a Mott 
insulator. 

*When the number of particles per period of the po- 
tential is not rational, there can be no commensuration 
and the system is again a conductor. The macroscopic 
behavior of such weakly-incommensurate systems is the 
subject of this paper. Excluding special cases (such as the 
"filled band" of one particle per potential well), the Mott 
state is degenerate. When the number of particles per 
potential well is = p/ 1, where p and q are co-primes, 
there are q different registries of the crystal relative to 



the potential. An important configuration is the soliton, 
which is the lowest-energy way of joining two domains of 
different registry. Although there are q-1 such combina- 
tions, we will assume that energetic considerations select 
a preferred sequence of neighboring domains, so that the 
possibilities reduce to " soliton" and " antisoliton" . Then 
for V close to p/q the structure may be regarded as a 
low-density gas of solitons. This is a weakly interacting 
system, even though the original particles are dense and 
their interactions are not small. 

We note that in the path integral representation the 
gas of quantum solitons is equivalent to a gas of classical 
lines in two dimensions, and thus describes the behavior 
of magnetic flux lines in a large area Josephson junction 
(i.e. a thin normal interface between two superconduc- 
tors), or steps on the surface of a crystal. 

Moving a soliton through the system shifts the Mott 
crystal by one potential well, and thus at z/ — p/q moving 
q solitons shifts p particles past a given point, returning 
the crystal to its original registry: a soliton is a fractional 
particle Q. 

Solitons repel each other with a short-ranged force; 
this can be approximated by treating the solitons as 
a gas of noninteracting free fermions, with the chemical 
potential 



(1) 



where e is the energy cost to introduce a soliton, n is 
the density of solitons, and m is the soliton mass. The 
last term is a quantum effect coming from the overlap of 
soliton wavefunctions, and is more important than the 
classical effect coming from the overlap of soliton strain 
fields when the intersoliton distance exceeds the soli- 
ton width ^. 

In equilibrium the chemical potential must be constant 
along the system, i.e. e-|-7r^?i^n^/2m = const. Depending 
on the magnitude of the constant this equation can have 
a solution with n real (the uniform soliton conductor) , or 
may have no real solution at all. The latter means there 
are no solitons in the system (Mott state, n = 0) and the 
condition of equilibrium is irrelevant as there is nothing 
to equilibrate. The Mott crystal plays a role of a vacuum 
of solitons/antisolitons. 

When the system is placed in an external field which 



provides a potential energy U(x) for the solitons, the con- 
dition of equihbrium becomes ||] 

fi + U{x) e + 7r^?i^n^/2m + U{x) = const (2) 

In contrast to the translationaUy-invariant case U — 
const, this equation allows solutions in which a region 
occupied by solitons can coexist with the soliton-free vac- 
uum. 

For example, let us look at a system of finite length 
Q < X < L placed in a field that exerts a constant force 
F on each soliton. If the original particles are charged, 
then solitons will be charged as well so that the prac- 
tical realization of such field would be the electric field. 
Equilibrium in such a field is only possible if the system 
is disconnected from external leads which we will assume 
to be the case. Then U{x) = —Fx, and (2) has a solu- 
tion of the form 71^ = {2mF/Tr^h'^){x - x^). For F > 
there is a soliton vacuum for x < Xg, and the solitons 
will occupy the x > Xg half-space with a sharp boundary 
between the two regimes located at x = Xs- However, 
solitons and antisolitons are oppositely charged, so that 
the same constant force will push them in opposite direc- 
tions, and there can be a second region for x < Xa where 
there is a finite density of antisolitons. We will associate 
n > with a density of solitons (an excess of the original 
particles) and n < with antisolitons (a particle deficit). 
Then the equilibrium distribution is 

{\ 2mF{x — Xs) lir'^'h^ x > Xs 

I ^ 1 
— \/2mF{Xa — x)/TT'^h X < Xa 

The values of x^ and Xa are constrained by the condition 
of conservation of particles, which is also a constraint on 
the integral of n(x). 

(i) Assume that for F the ground state is the Mott 
crystal with no solitons present. A nonzero field F will 
promote the creation of soliton- antisoliton pairs. The 
change of the system energy upon creation of a soliton- 
antisoliton pair is given by e{y) = 2e — Fy where the first 
term is the energy cost to create two solitons of oppo- 
site kind separated by a distance y ^ £, while the second 
term is the energy gain in external field. The pairs for 
which e(y) < 0, i.e. those separated by a distance bigger 
than 2e/F will be present in equilibrium. As a result 
the external field F polarizes the Mott insulator by cre- 
ating and spatially separating solitons and antisolitons. 
If the field F pushes solitons to the right and antisoli- 
tons to the left, the density distribution n(x) is odd in x 
about a; = 0, and we have n = for < x < e/F and 
n = {2mF/TT^h^y/^{x - e/F)^/^ for x > e/F : a strip 
of Mott phase of width 2e/F separates solitons from an- 
tisolitons. This conclusion is only true for a sufhciently 
large system whose size L exceeds the size of the Mott 
strip 2e/F; otherwise pair creation is not profitable and 
the vacuum is the lowest energy state. As the field in- 
creases, the Mott strip shrinks, and our description fails 



at very large fields of order e/^ when the size of soliton- 
free region becomes comparable with the soliton width ^. 

(ii) Assume that for F = the ground state is a soliton 
conductor. A sufficiently small nonzero field F pushing 
solitons to the right will turn the uniform soliton dis- 
tribution into n = (2TOF/7r^7i^)^/^(a; — Xs)^^"^ illustrated 
in Figure la (for this case Xs is outside of the physi- 
cal region). As the field increases, at some F there will 
be a marginal configuration (xs — 0) (shown in Figure 
lb) where the soliton density vanishes at the left end 
of the system. At larger fields {xs > 0) a soliton-free 
Mott strip forms at the left end of the system (Figure 
Ic) . Upon further increase of the field when solitons get 
pushed sufficiently far away from the left, there will be 
another marginal configuration for which the size of the 
Mott strip is exactly 2e/F. At largest fields, antisolitons 
nucleate at the left end of the system {xa > 0), the num- 
ber of solitons at the right increase by the same amount, 
and the size of the Mott strip stays equal to 2e/F there- 
after (Figure Id). 
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Soliton distribution in a soliton conductor for various 
applied fields, (a) For small fields, the soliton distribu- 
tion is nonuniform but everywhere nonzero, (b) There is 
a critical value of the field, where the soliton density is 
zero at one end of the system; (c) for larger fields, there 
is a region where n{x) =0. (d) For even larger fields, 
soliton-antisoliton pairs nucleate. 

The equilibrium configurations of solitons in an exter- 
nal field are related to the equilibrium crystal shapes of 
three- dimensional crystals, since the latter system can 
be viewed as equilibrium of atomic steps 1^. 

We now turn to a discussion of nonequilibrium effects, 
and hereafter we assume that there arc no external fields 
present. The force exerted on a given soliton by its neigh- 
bors is —d^i/dx, which will cause it to drift with the ve- 
locity V = —^dpujdx, where 7 is phenomcnological fric- 
tion constant. The resulting current of solitons is 



j = nv = —•yndfi/dx = —Tdj-i/dx = —bn^dn/dx (4) 

where F — "/n is the system mobihty, b = jn'^h'^ /m is 
a dynamical parameter, and we used Eq.(l) to compute 
(9/x/9x. The mobihty F is hnear in the sohton density n, 
which imphes that the conductivity of the system van- 
ishes hnearly in deviation from commensuration as the 
Mott insulator phase is approached from the conductor 
side. This transport theory explains the linear drops in 
resistance seen by van Oudenaarden et al. in the vicin- 
ity of the Mott phases, and parallels the flux line mech- 
anism of resistivity of a type-II superconductor in the 
vortex state pO| , and the growth regime of vicinal crys- 
tal surfaces via the motion of steps |^ . 

Conservation of soliton number within a region implies 
a continuity equation, which provides an equation of mo- 
tion for n(x, t): 

dn/dt^~dj/dx^b-^{n^^) (5) 
ox ox 

When the function n{x,t) changes sign (both solitons 
and antisolitons present) the relaxation within the two 
regions continues to be described by Eq.(5); however, the 
behavior at the interface n = where the solitons and 
antisolitons annihilate requires further discussion. 

When a soliton- antisoliton pair annihilates, an energy 
2e is released [first term of (1)]. Viewed classically, the 
force that corresponds to this potential energy has only 
the range of order ^, since solitons do not communicate 
beyond this scale. From our macroscopic viewpoint this 
would imply no interaction between solitons of opposite 
kind at all, and the dynamics of annihilation would be de- 
termined by the rate at which they get pushed together 
by the excess of particles elsewhere. This overlooks a 
quantum effect: there always is a finite annihilation prob- 
ability of a pair whose size is bigger than ^. This possibil- 
ity of future annihilations gives an effective longer- range 
attraction between solitons of opposite type. However, 
the dominant effect associated with the region where soli- 
tons and antisolitons annihilate is that the mobility van- 
ishes [Eq.(4)], so that attractive forces between solitons 
and antisolitons have no important macroscopic effect on 
dynamics, and (5) is still valid. 

Equation (5) resembles the diffusion equation, but the 
"diffusion constant"!) = bn^ is density-dependent. This 
equation has previously appeared in the theories of shock 
waves, filtration ||Tl[| , and dynamics of crystal surfaces 
@- 

The most remarkable property of (5) is that it allows 
self- similar moving solutions with a sharp front where 
the soliton density vanishes. Indeed, assume there is such 
a boundary. During a sufficiently small time interval the 
velocity of the front v can be taken approximately con- 
stant. The distribution of soliton density near the front 
can be sought in the form n(x — vt). Substituting this 



into (5) and integrating twice between the front position 
Xf[n{xf) — 0] and x we find 

vn'^/2-jfn + (jf/v) \n\{vn + jf) / jf \ = v'^{xf - x)/b 

(6) 

where j / has a meaning of the current density at the front 
position X — Xf. The shape of the density profile near 
the front edge depends on whether jf is zero or finite: 

(i) The case jf = has been considered previously . 
Taking in (6) the limit jf = we find = (2v/b){x f — x) 
which implies that the solitons are present in the region of 
space satisfying v{xf — x) > 0. The region v{xf — x) < 
is soliton- free which is consistent with the assumption 
that the current density jf vanishes at the front position. 
The density distribution near the front is given by 

n = ±y^fc3 (7) 

This type of front describes a cloud of soli- 
tons/antisolitons invading a Mott phase. As in the equi- 
librium case (3), the front has a square-root singularity 
(7). The steepness of the moving front is determined by 
the velocity v of the front. 

(ii) If n(xf) = is the annihilation front then solitons 
and antisolitons get pushed by their neighbors towards 
X = Xf, and the current density jf must stay finite there. 
As far as we can tell, this type of shock has not been stud- 
ied previously. Solving (6) to lowest nonvanishing order 
in \vn/ jf \ <C 1 the density distribution near the front can 
be found as 

n = [3jf{xf - x)/b]'/^{l + iv/Ajf)[3jfixf - x)/b]'^'} 

(8) 

This type of front describes solitons invading an antisoli- 
ton conductor (or the other way around). The peculiar 
feature of (8) is that the magnitude of the force acting on 
a soliton dfx/dx cx ndn/dx cx {xf — x)~^^^ oc n~^diverges 
at the annihilation edge. This divergence offsets the van- 
ishing of the system mobility and thus overcomes a dy- 
namical bottleneck. As a result the current density (4) 
stays finite at the annihilation edge which is the ultimate 
reason for the leading cubic-root singularity in (8). 

As Eq.(8) shows, there can be both moving and sta- 
tionary annihilation fronts. In fact, for v = 0, Eq.(8) 
provides the exact steady-state solution to (5). The pa- 
rameter jf will now have the meaning of a fixed current 
fiowing through the system which is the same everywhere. 

The front coordinate x/, velocity v = dxf/dt, and 
annihilation current jf entering Eqs.(6)-(8) are undeter- 
mined functions of time which should be found from the 
complete solution of the problem. For the special case 
of a constant velocity front, we have Xf = vt+ const, 
and (7) becomes an exact solution to the problem of a 
constant velocity soliton invasion of the Mott phase. Sim- 
ilarly, Eq.(6) provides an exact solution to the problem 



of a constant velocity annihilation front which can be 
visualized as follows: 

Assume v > so that the front moves to the right. 
Far to the left of the front position, x <C Xf, the 
front profile is approximately given by (7). Far to the 
right, X ^ Xf, the density tends to a constant value 
rioo = ~iflv, E^nd then the solution (6) can be approx- 
imated by n = noo{l — exp[(t)/6n^)(a;/ — x)]}. This is 
exactly the density profile one would get ahead of a con- 
stant velocity front if the relaxation could be described 
by the linearized version of (5) with the diffusion con- 
stant D = hri^; the size of the perturbed region ahead of 
the front is given by the diffusion length D/v — hri^/v. 

Below we consider a few more practically relevant ex- 
amples when the front velocities are not necessarily con- 
stants. 

(i) Imagine a Mott insulator for which extra particles 
have been added to a region of the system. This can be 
described as a segment of nonzero soliton density imbed- 
ded between two half- infinite commensurate domains. 
The number of particles added determines the total num- 
ber of solitons N, which is conserved: n(x, t)dx = N. 
Arbitary initial particle distributions will relax asymptot- 
ically to a common form 

r.(x, t) = ( Jvl ) y^i_^^.2/4^V6^ (g) 

which has the characteristic square-root singularity (7) 
at its edges. 

(ii) Consider the case that the average density of par- 
ticles is appropriate for commensuration but the initial 
distribution is inhomogeneous. For example, we could 
have a Mott insulator everywhere except in a small re- 
gion, where there is an excess of particles on one side, 
and a deficit on the other. In the soliton description, 
the two regions will contain solitons of different types 
-solitons and antisolitons - which can be described by 
a single function n(x). Here are two examples in which 
soliton-antisoliton annihilation plays a role: 

Consider a periodic distribution of solitons of both 
signs such that the total soliton "charge" is zero. The 
spatial periodicity of the distribution will be preserved by 
time evolution; this implies that the asymptotic solution 
to Eq.(5) should be sought in the form |l2| n = t^°'f{x). 
It is readily seen that the equation of motion (5) deter- 
mines the exponent to be a = 1/2. The density decays 
because solitons succeed in getting through the zero mo- 
bility region, by having a singularity in the applied force: 
the function f(x) has Ax^^^ singularities at every point 
where it changes sign [ p^ in agreement with the general 
argument leading to (8). 

A related example is a soliton distribution such that 
n{—x) = —n{x). This symmetry will be preserved by 
the time evolution. Due to annihilation events at a; = 0, 



the total number of solitons of either type is not con- 
served; however, the dipole moment P = xn{x, t)dx 
is conserved Q by the equation of motion (5). The ex- 
pansion rate is smaller than in previous case of the spread 
of a region of excess solitons because of annihilation of 
soliton-antisoliton pairs at x = 0; the envelope decays 
more slowly than in the periodic case because the soli- 
tons are not confined. The analytic solution for this case 
is 

-M-^{^f"[^-{fjt"f'' (10) 

where A and B are constants. At the extremal edges, 
there is again the V Ax singularity, (7), but near the ori- 
gin •n{x) cx x^/'^ [see Eq.(8)], as is required to have a net 
current of solitons through the mobility bottleneck. 
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